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Summary 

The NASA Storm Hazards Program was initi- 
ated to improve the knowledge and understanding 
of atmospheric processes and their effects on aircraft 
flight operations in the vicinity of storms. An instru- 
mented F-106B airplane was used to penetrate thun- 
derstorms and to make measurements of the wind at 
points along the flight path. During these flights, a 
ground-based Doppler radar also made independent 
measurements of the wind. This report presents the 
results of five airplane penetrations of two storms in 
1980 and six penetrations of one storm in 1981. Com- 
parisons were made between the wind velocity mea- 
sured by the radar and the airplane. Least-squares 
regression lines were calculated, and correlation co- 
efficients of 0.88 and 0.78 were obtained for the 1980 
data and part of the 1981 data, respectively. This 
combined data set also had a correlation coefficient 
of 0.88. Although the correlations may be as good 
as can be expected for comparisons of this type, it 
is suggested that improvements in the experimental 
technique could possibly improve the results. 

Introduction 

The Storm Hazards Program was initiated to im- 
prove the knowledge and understanding of atmo- 
spheric processes and their effects on aircraft flight 
operations in the vicinity of severe storms. Results 
from the program will be useful in improving the de- 
sign of aircraft and systems for storm hazards pro- 
tection and for the detection and avoidance of haz- 
ardous weather conditions (ref. 1). Since its in- 
ception, the primary emphasis of the program has 
been on characterizing lightning and its effects on 
advanced aircraft avionics and composite structures 
(refs. 2 and 3). In the initial phase of the program, 
an instrumented NASA twin-engine light transport 
airplane flew around the periphery of storms to de- 
tect and locate lightning discharges. Later, thunder- 
storm penetrations were made in an F-106B airplane 
instrumented to make extensive measurements of the 
thunderstorm environment. These measurements in- 
cluded the three wind components at points along the 
flight path. Also during these flights, a ground-based 
Doppler radar made independent measurements of 
the wind velocity along the radar beam. Interest in 
using a Doppler radar for characterizing the wind ve- 
locity in the vicinity of thunderstorms has increased 
recently because of its potential application for im- 
proving aviation safety. For this reason, these two 
independent measurements of the wind velocity in 
thunderstorms were compared. The purpose of this 
paper is to present the results of those comparisons. 


Flight operations in 1980 and 1981 were con- 
ducted in Norman, Oklahoma, in a cooperative effort 
with the National Severe Storms Laboratory (NSSL). 
Generally, the method of operation was to locate a 
storm by means of ground-based radars, deploy the 
NASA F-106B airplane, direct the airplane to the 
vicinity of the storm, and fly through the storm sev- 
eral times before returning to the base. No attempt 
was made, however, to achieve simultaneous or co- 
extensive velocity measurements by the airplane and 
ground-based Doppler radar. Flight operations in- 
cluded five thunderstorm penetrations in June 1980 
and six penetrations in April 1981. Radial velocity 
contours measured by Doppler radar are presented 
for each thunderstorm penetration with the F-106B 
ground track superimposed on the contours. Com- 
parison plots of the airplane and radar radial velocity 
data are presented for 1980, 1981, and the combined 
data for both years. 

Symbols and Abbreviations 

A z vertical acceleration, positive up, 

g units 

AIS aircraft instrument system 

a unit vector from airplane to radar 

g acceleration due to gravity, 9.807 m/s 2 

h p altitude determined from pressure 

measurement (pressure altitude), m 

INS inertial navigation system 

l distance between INS package and 

flow vanes along longitudinal axis of 
F-106B, m 

NSSL National Severe Storms Laboratory 

p roll velocity of F-106B about body 

x-axis, rad/s 

q pitch velocity of F-106B about body 

y-axis, rad/s 

q f c measured dynamic pressure 

r yaw velocity of F-106B about body 

z-axis, rad/s 

T integration interval, s 

t time, s 

V velocity of F-106B relative to air mass, 

m/s 

V i inertial velocity of F-106B, m/s 

V X ’) V y ' } V Z ' components of V/, m/s 


( Vz')o 

initial value of inertial velocity used in 
equation (A6), m/s 

W 

wind velocity, m/s 

W r 

radial wind velocity, m/s 

Wra 

radial wind velocity computed from 
F-106B data, m/s 

W rr 

radial wind velocity computed from 
Doppler radar data, m/s 

w x .,w y ., 

W z > components of W in inertial axis 
system, positive north, east, and 
down, m/s 

x,y,z 

body axes of F-106B, dimensionless 

x',y',z' 

inertial axes, dimensionless 

a 

angle of attack as measured by flow 
vanes and corrected for flow-field 
effects, deg 


angle of attack as measured by flow 
vanes, deg 

P 

angle of sideslip as measured by flow 
vanes (no flow-field corrections used), 
deg 

r 

transformation matrix defined by 
equation (A3) of appendix 

A p s 

difference between measured static 
pressure and free-stream static 
pressure 


Euler angles, deg 


Experimental Apparatus 

The experimental apparatus used in the Storm 
Hazards Program at the NSSL is described in detail 
in references 1-5. For that reason, only the equip- 
ment and systems pertinent to data presented herein 
are described in this paper. 

Research Airplane and Systems 

A photograph of the NASA F-106B airplane is 
shown in figure 1. The airplane is an all-weather, 
two-seat, delta-wing interceptor with a leading-edge 
sweep of 60°. Modifications were made to the air- 
plane especially for the Storm Hazards Program and 
were described in reference 3. Basic characteristics 
of the airplane (from ref. 2) are listed in table I. The 
airplane was equipped with an X-band color weather 
radar to display precipitation reflectivity contours, 
an aircraft instrument system (AIS), and an inertial 
navigation system (INS). 


The aircraft instrument system measurements 
included angles of attack and sideslip, airplane 
body rates, pitot-static and dynamic pressures, to- 
tal temperature, airplane center-of-gravity accelera- 
tions, and lateral and longitudinal stick position. A 
14-track magnetic tape recorder and time code gen- 
erator were used for recording time histories of the 
various measurements. The AIS also recorded the 
output of the INS. Data were recorded at a rate of 
20 samples per second. 

The INS was designed for commercial use to pro- 
vide aircraft navigation. The INS measurements in- 
cluded latitude and longitude; pitch, roll, and true 
heading; vertical acceleration; and horizontal com- 
ponents of airplane ground speed. Experience with 
this system as used on the F-106B missions indi- 
cates a typical position error of about 1 nautical 
mile per hour of flight. 

A list of the parameters recorded by the AIS and 
used to derive the three components of wind are pre- 
sented in table II. Also listed are the instrument 
range, frequency response, lcr error, and measure- 
ment resolution. The dynamic-pressure sensor and 
flow vanes were located on the nose boom. Rate gy- 
ros were located in the internal weapons bay. The 
INS and the longitudinal, lateral, and normal ac- 
celerometers were mounted at the center of gravity 
of the airplane. The vertical accelerometer was an 
integral part of the INS platform. Pitot-static pres- 
sure and air total temperature probes were located 
beneath the nose of the airplane. 

Doppler Radar 

The NSSL has a matched pair of Doppler radars, 
one located at Norman and the other at Cimarron, 
for scanning thunderstorms. During the Storm Haz- 
ards Program, no storms entered the coverage area 
common to both radars, and for that reason, all 
Doppler radar data were obtained by the radar lo- 
cated at Norman. 

Doppler radar characteristics for the radar lo- 
cated at Norman are listed in table III. (See ref. 6.) 
The radar had a half-power beam width of 0.8°, a 
wavelength of 10.5 cm, and a peak power of 750 kW. 
The pulse width was 1 /xs, and a scan rate of 6 deg/s 
was used. The pulse repetition time was normally 
768 /xs, and 64 samples were used to estimate the 
mean velocity. The dwell time was typically 50 ms. 
After completing an azimuthal scan, the elevation 
was stepped and the scan was repeated. A scan at 
each of three elevation angles was made, and data 
at three altitudes were obtained such that the inter- 
val between the low and high altitudes included the 
altitude at which the F-106B was flying. 
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The volume in space sampled by the radar (radar 
resolution volume) for a zero antenna scan rate has 
dimensions defined by the pulse width (depth of 
volume), antenna beam width, and the range from 
the radar to the sample volume. (Beam width and 
range define the cross-sectional area perpendicular to 
the centerline of the beam.) The resolution volume 
is an inherent characteristic of the radar and is 
smaller than the actual volume sampled in space. 
The reason for this is that the Doppler spectrum is 
determined from many transmitted pulses, and as the 
radar scans, the volume contributing to the spectrum 
is larger than the resolution volume. This larger 
volume is defined to be the measurement volume 
(ref. 7). The measurement and resolution volumes 
were assumed to be equal for these experiments. 

Analysis Method 

F-106B Measurements 

The wind radial velocity along the Doppler radar 
beam from the F-106B to the radar was derived from 
measurements of the parameters listed in table II. 
The wind vector was determined by the method 
presented in the appendix. The wind velocity data 
as determined from the F-106B measurements were 
computed at a rate of 20 samples per second. Since 
the Doppler radar radial velocity was actually an 
average over the radar measurement volume, the 
airplane data were averaged over a 5-s interval in 
order to make the comparisons between the radar and 
airplane data. This interval was the time required 
for the F-106B to travel about 1 km, which is the 
diameter of a sphere with a volume equivalent to 
the radar measurement volume at the nominal radar 
range (140 km) of the tests. The sensitivity of the 
airplane radial velocity measurements to averaging 
intervals of 2, 5, and 10 s is shown in figure 2. 
As would be expected, larger averaging intervals 
smooth the data by removing variations in the radial 
velocity due to turbulence and sharp-edged updraft 
or downdraft regions. 

Doppler Radar 

Doppler radar measures target speeds of precip- 
itation particles as the beam sweeps the thunder- 
storm. Motion of the particles toward or away from 
the radar (radial velocity) is deduced from minute 
shifts in the microwave frequency of precipitation 
echoes. Because the volume sampled in space by 
the radar contains many targets, the radar actually 
measures a spectrum of radial velocities. The first 
moment of the Doppler spectrum is the mean radial 


velocity of the particles in the volume sampled. Prin- 
ciples of pulsed Doppler radar may be found in ref- 
erences 8 and 9. Some details on Doppler radar data 
reduction procedures may be found in reference 1. 

Doppler radar radial velocity contours on con- 
stant altitude surfaces as a function of position from 
the Doppler radar as shown in figure 3 were pro- 
vided by the NSSL. The contours were produced 
from sector-scans of the storm that were made dur- 
ing the time that the F-106B was flying through the 
storm. Generally the radar data measurements re- 
quired about 3 min and were synchronized with the 
time of the F-106B measurements within 2 min. Ve- 
locity contours for the 1980 data were provided by 
the NSSL at the altitude flown by the F-106B. Ve- 
locity contours for the 1981 data, however, were pro- 
vided at three altitudes, and it was necessary to in- 
terpolate to the F-106B altitude. The interpolation 
method used is illustrated in figure 3. 

Velocity contours at three altitudes are shown in 
figure 3(a). It was assumed that the wind fields were 
frozen during the time required for the F-106B to 
traverse the thunderstorm. There are several pos- 
sible approaches for comparing the Doppler frozen 
wind field data, which were measured as a function 
of position, with the airplane data, which were mea- 
sured as a function of time along the flight path. One 
possible approach, for example, would be to super- 
impose the velocities measured by the airplane at 
positions on the frozen wind fields. A second ap- 
proach would be to take the Doppler radar contour 
data at a given position and transform it to an ef- 
fective time history along the flight path. The lat- 
ter approach was chosen to make these comparisons. 
For this method a radial velocity effective time his- 
tory was constructed at each altitude to represent 
the Doppler radar radial velocity at points along the 
F-106B flight path. This was done by calculating the 
F-106B ground track as a function of position from 
the Doppler radar and plotting the track over the 
contours at each altitude. Elapsed time along the 
flight path is indicated in figure 3(a). The radial ve- 
locity values as a function of position were extracted 
at each altitude and plotted in figure 3(b) as a func- 
tion of elapsed time along the flight path. For the 
penetration shown, the F-106B altitude was 3.7 km. 
Values of the radial velocity at that altitude were de- 
termined by interpolation and are also shown in fig- 
ure 3(b). Direction of flight through the storm and 
direction to the Doppler radar at Norman from the 
F-106B are also indicated in figure 3(a). 

Results and Discussion 

Radial velocity data derived from Doppler radar 
and F-106B measurements for flights made in June 
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1980 are shown in figures 4-8. Figures 4 and 5 show 
results for two passes through the thunderstorm at 
an altitude of 4.6 km on June 16. Figures 6-8 show 
results for three passes through another thunder- 
storm at an altitude of 4.9 km on June 17. The 
average distance from the radar to the F-106B for 
the data shown in figure 4 was 113 km, and for 
the data shown in figure 7, the average distance was 
155 km. Average distances for the other penetrations 
fell within those boundaries. 

Ground tracks derived from the F-106B INS data 
are shown in the (a) parts of figures 4-8 superim- 
posed on the Doppler radar radial velocity contour 
plots. The direction of flight of the F-106B airplane 
through the storm is indicated on each figure as well 
as the direction toward the Doppler radar at Norman. 

The (b) parts of figures 4-8 show time histories 
of the radial velocity derived from the F-106B mea- 
surements (see eqs. (A2) and (A4) in the appendix) 
and from the Doppler radar data presented in the (a) 
parts of figures 4-8. These data were used to make 
the comparison plot shown in figure 9. The regres- 
sion line in figure 9 was derived from 96 data points 
obtained by reading the radial velocity values at 15-s 
intervals along the ground tracks in figures 4-8 and 
matching these points in time with the F-106B val- 
ues. The correlation coefficient was 0.88. 

Figures 10-15 show the radial velocity contour 
plots and time histories for a flight conducted on 
April 22, 1981. The F-106B made six passes through 
the thunderstorm: two passes at an altitude of 

3.7 km, one pass at 4.6 km, two passes at 6.1 km, 
and one pass at 6.0 km. Contour plots at the 
altitude nearest the F-106B altitude are shown in 
figures 10-15. These altitudes were within 400 m of 
the F-106B altitude for the data in figures 10-14, 
and within 300 m for the data in figure 15. The 
range from the radar to the F-106B generally was 
greater than 120 km. Figure 16 shows a comparison 
plot of the 1981 radial velocity data derived from 
figures 10-15. The regression line was derived from 
93 data points, and the correlation coefficient was 
0.47. 

Comparisons of the 1980 Doppler radar and air- 
plane radial velocity data summarized in figure 9 
were considerably better (Correlation coefficient = 
0.88) than the 1981 data summarized in figure 16 
(Correlation coefficient = 0.47). It is reasonable to 
expect that the correlation coefficients would be com- 
parable. For example, the 1980 data were derived 
from two flights (five passes) through two thunder- 
storms, and the 1981 data were derived from one 
flight (six passes) through one thunderstorm. Meth- 
ods used to derive the Doppler radar radial velocity 
time histories along the F-106B flight path for both 


data sets were essentially the same. Also, both re- 
gression lines were calculated from about 100 data 
points. 

Apparently the correlation coefficient in figure 16 
was lower than the correlation coefficient in figure 9 
because of significant differences in the F-106B and 
Doppler radar radial velocity values for three of the 
six passes (i.e., passes 3, 4, and 6) through the thun- 
derstorm in April 1981. These differences were shown 
in figures 12(b), 13(b), and 15(b). It is possible that 
the differences were caused by large variations in the 
wind field within the radar measurement volume. For 
example, the radar measurements were averages over 
the measurement volume, whereas the F-106B mea- 
surements were along the flight path (a line) in a hor- 
izontal plane in the measurement volume. The aver- 
age distance of the F-106B from the Doppler radar 
for each of the three passes was about 140 km. This 
means that the cross section of the radar measure- 
ment volume perpendicular to the beam had a diam- 
eter of about 2 km, and the radar would have aver- 
aged the wind field with that cross section. It is also 
possible that the radar-derived altitudes for the con- 
tour plots could at times have been in error at these 
long ranges because of refractivity variations along 
the propagation path. In all cases, the radar data 
were subjected only to a standard earth curvature 
refractivity correction. 

Figures 17(a)-17(c) show the F-106B radial veloc- 
ity time histories for penetrations 3, 4, and 6 derived 
by interpolating between three altitudes. These fig- 
ures show that the time histories above and below 
the F-106B altitude were similar, and radial veloc- 
ity magnitudes were not significantly different, ex- 
cept for the first half of penetration 4 in figure 17(b) 
and the latter part of penetration 6 in figure 17(c). 
To make these interpolations, it was assumed that 
the airplane flew through the storm at a constant 
altitude that was equivalent to the altitude at en- 
try to the storm. Actually the altitude within the 
storm may have changed because of sudden changes 
in the wind field. Experience has shown that these 
changes rarely exceed 200 m during a penetration but 
may have occasionally been 300 m (approximately 
1 time per 100 penetrations). Also the altitude at en- 
try was based on barometric pressure and may have 
been different from the geometric altitude represent- 
ing the radar data at that altitude. These uncer- 
tainties along with probable time differences between 
the airplane and radar measurements, which could 
mean that the storm dynamics at a given altitude 
had changed, may account for the trends shown in 
figures 17(b) and 17(c) and for the differences seen in 
figures 13(b) and 15(b). These trends (see figs. 17(b) 
and 17(c)) suggest that the differences in the wind 
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field as a function of altitude may account for the 
smaller correlation coefficient for the 1981 data in 
figure 16. In fact, when the data for penetrations 3, 
4, and 6 were deleted from the 1981 data set, the 
results were as shown in figure 18. The comparisons 
improved, and the correlation coefficient increased 
from 0.47 to 0.78. 

Figure 19 shows comparison plots for the 1980 and 
1981 data combined. Data presented in figure 19(a) 
include all the thunderstorm penetrations, whereas 
the data presented in figure 19(b) exclude penetra- 
tions 3, 4, and 6 of the 1981 data. As might be 
expected, the comparisons were improved by exclud- 
ing the three penetrations, and the correlation coef- 
ficient was increased from 0.78 to 0.88. Figure 19(b) 
may represent the best correlation that can be ex- 
pected between Doppler radar and airplane data for 
similar experimental conditions as presented here. It 
is suggested, however, that several improvements in 
the experimental technique may improve the results. 
These improvements might include, for example, ac- 
counting for airplane position errors and making the 
airplane and radar measurements at the same time 
and close to the same position in space by slaving the 
radar to the airplane position. These improvements 
would properly locate the airplane position within 


the storm and minimize differences due to thunder- 
storm dynamics. 

Concluding Remarks 

A comparison of wind velocity measurements 
made in thunderstorms by an instrumented airplane 
and a ground-based Doppler radar was made. Com- 
parisons were made for five airplane penetrations 
of two thunderstorms in 1980 and six penetrations 
of one thunderstorm in 1981. No attempt was 
made to achieve simultaneous or coextensive mea- 
surements, and consequently, an interpolation was 
necessary to match the two measurements in space 
and time. Correlation coefficients for the 1980 data 
and part of the 1981 data were 0.88 and 0.78, re- 
spectively. The correlation coefficient for this com- 
bined data set was 0.88. Some of the 1981 data 
showed large discrepancies between the airplane- and 
radar-measured data. Potential contributing fac- 
tors for the discrepancies were identified and dis- 
cussed. These correlations may represent the best 
that can be expected for this type of comparison, but 
it is suggested that improvements in the experimen- 
tal technique would probably improve the results. 

NASA Langley Research Center 
Hampton, VA 23665 
February 1, 1985 
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TABLE I. CHARACTERISTICS OF F-106B RESEARCH AIRPLANE 


Length, m (ft) 

Height, m (ft) 

Wing span, m (ft) 

Wing area (gross), m 2 (ft 2 ) . . . 

Wing chord at root, m (ft) . . . 

Aspect ratio 

Wing sweepback angle 

Empty weight, N (lbf) 

Gross take-off weight, N (lbf) . . . 

Engine 

Thrust (military) at sea level, N (lbf) 
Maximum thrust, N (lbf) 


21.54 (70.67) 

6.17 (20.25) 

11.70 (38.39) 

64.83 (697.82) 

9.07 (29.77) 

2.118 

60°6'13" 

.... 116 543 (26 200) 

.... 160 710 (36 129) 

J75-P-17 axial flow turbojet 
.... 71 616 (16 100) 

.... 108 981 (24 500) 


TABLE II. WIND AND TURBULENCE MEASUREMENTS ON F-106B AIRCRAFT INSTRUMENTATION SYSTEM 


Measurement 

Instrument 

range 

l cr error 

Resolution 

Required flat 
to frequency 
response, Hz 

Notes 

Static pressure, psia* 

0 to 15 

0.07 

0.0004 

3 

Quartz transducer — under 
nose Pitot-static head 

Dynamic pressure, psi* .... 

0 to 3 

0.03 

0.0004 

6 

Quartz transducer — nose 
boom for total pressure 

Angle of attack, deg 

±30 

0.09 

0.015 

10 

Flow vane on nose boom 

Angle of sideslip, deg 

±30 

0.09 

0.015 

10 

Flow vane on nose boom 

Body rates, rad/s 

±1 

0.01 

0.005 

5 

Mounted in weapons bay 

Longitudinal and lateral 

acceleration, g units .... 

±1 

0.01 

0.004 

10 

Mounted at center of gravity 

Normal acceleration, g units . . 

6 to -2 

0.01 

0.012 

10 

Mounted at center of gravity 

Total air temperature, °C . . . 

35 to —75 

0.03 

0.205 

0.5 

Probe beneath nose 

Pitch attitude, deg 

±30 

0.09 

0.03 

5 

INS platform — continuous 
synchronization 

Roll attitude, deg 

±80 

0.24 

0.08 

5 

INS platform — continuous 
synchronization 

Ground speed, m/s 

±1248 

0.9 

0.0003 

3 

INS platform — digital 

Latitude 

±90° 

1/60 deg/hr 

0.14" 

0.2 

INS platform — digital 

Longitude 

±180° 

1/60 deg/hr 

0.14" 

0.2 

INS platform — digital 

TVue heading 

0° to 360° 


0.14" 

3 

INS platform — digital 

Vertical acceleration, g units . . 

3 to -1 

0.005 

0.008 

5 

INS platform acceleration 


1 psi = 6.9 kPa. 









TABLE III. NSSL DOPPLER RADAR CHARACTERISTICS 


Antenna: 

Shape 

Diameter, m 

Half-power beam width, deg 

Gain, dB 

First side lobe level, dB . . 

Polarization 

Root-mean-square surface deviation, mm . . . 

Transmitter: 

Wavelength, cm 

Frequency, MHz 

Peak power (typical), kW 

Pulse width, /z s (m) 

Pulse repetition time (typical), /zs 

Receiver: 

System noise (typical),* dB 

Transfer function 

Dynamic range (typical), dB 

Band width, MHz 

3 dB 

6 dB 

Noise power (typical), dBm 

Doppler time series data acquisition: 

No. of simultaneous range gates 

No. of blocks along radial 

Range gate spacing, m 

Azimuthal sample spacing, deg 

Automatic elevation increment, deg 

No. of samples in time series 

Complex video digital word length 

Mean velocity and spectrum width data acquisition: 

No. of range gates 

Range gate spacing, m 

No. of samples in estimate 

Recorded word length 

Velocity 

Width 

Intensity data acquisition: 

No. of range gates 

Range gate spacing, m 

No. of samples in estimate 

Time window function 

Recorded word length 

*Measured at receiver calibration port. 


Parabolic 
. 9.15 

. 0.81 
. 46.8 

. -21 
Vertical 
2.8 


10.52 

2850 

750 

1 (150) 

768 

4 

. . Doppler — linear 
Intensity — logarithmic 
80 

0.60 

0.85 

-108 


16 

8 

150, 300, 600, 1200 
. . . . 0.1 to 9.9 

. . . . 0.1 to 9.9 

• 2"; n = 1,2, ...,13 
12 bits (binary) 


.... 762 

. . . . 150 

2 n ;n = 5,6, 7, 8 


6 bits (2’s complement) 
6 bits, 0-32 (binary) 


762 

150,600 

. 2 n — 1; 2 n ; n = 3,4 
Linear or exponential 
. . . 6 bits (binary) 
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TABLE III. Concluded 


General features: 

Azimuthal antenna rotation rate, deg/s 0.1 to 10 

Maximum unambiguous range, km 115 

Maximum unambiguous velocity (typical), ms -1 ±34.2 

Plan-position indicator capability Yes 

Range-height indicator capability • Yes 

Coplane capability No 

Calibration and timing: 

Radome loss (two-way), dB 2.0 

Waveguide loss (two-way, from calibration port 

to antenna feed assembly), dB 2.0 

Transmitter/receiver circulator loss (receiver port), dB 1.7 

Receiver detection loss,* dB 2.7 

Signal generator cable loss, dB 9.5 

Directional coupler loss, dB 36.2 

Delay to transmitter pulse,* /is 5.4 

Pulse pair processor integrator delay, § /is 3.5 

Waveguide delay, /is 0.1 

Doppler sample-hold aperture, /is 0.4 

Doppler sample-hold time constant, /is 0.15 

Pulse pair processor integrator sample-hold aperture, /is 0.04 

Pulse pair processor integrator sample-hold time constant, /is 0.01 

Receiver initial delay, /is 1.3 

Receiver nominal delay, /is 1.8 

Receiver rise time, /is 1.0 


*See memo, “Reflectivity Equation for NSSL’s WDS-71 10 cm Doppler Radar,” Feb. 21, 1973, R. J. Doviak 
and D. Sirmans. 

*Delay from system pretrigger to first sample-hold operation. 

§Pulse pair processor integrator delay — expanded, 6.5 /is. 
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Figure 1 . NASA F- 106 B research airplane used in Storm Hazards Program 
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Figure 2. Typical radial velocity time histories derived from F-106B measurements with averaging intervals of 
2, 5, and 10 s. 
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( a ) F-106B ground track superimposed on Doppler radial velocity contours at three altitudes. Contours are 
labeled in meters per second. 

Flg ure 3. Derivation of Doppler radar radial wind velocity at points along F-106B flight path by superposition 
and by interpolation between horizontal planes containing Doppler-measured radial velocity contours. 







(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



(b) Radial velocity. 

Figure 4. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity at 
points along flight path. Penetration 4; June 1980. 
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Nominal altitude = 4.6 km 
F-106B altitude * 4.6 km 



(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



(b) Radial velocity. 

Figure 5. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity at 
points along flight path. Penetration 5; June 1980. 
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Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 


(b) Radial velocity. 

Figure 6. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity at 
points along flight path. Penetration 1; June 1980. 
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100 


Nominal altitude = 4.9 km 
F-106B altitude = 4.9 km 



Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 

(b) Radial velocity. 

Figure 7. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity at 
points along flight path. Penetration 2; June 1980. 
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no 


Nominal altitude = 4.9 km 
F-106B altitude = 4.9 km 



Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 


00 



a: 


Elapsed time, min 

(b) Radial velocity. 

Figure 8. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity at 
points along flight path. Penetration 3; June 1980. 
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O Penetration 4, figure 4 
□ Penetration 5, figure 5 
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Figure 9. Comparison of Doppler radar and F-106B radial velocity data from June 1980. 
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Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



(b) Radial velocity. 

Figure 10. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 1; April 1981. 
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(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



(b) Radial velocity. 

Figure 11. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 2; April 1981. 
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Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 

(b) Radial velocity. 

Figure 12. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 3; April 1981. 
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Distance east of radar, km 

(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 


(b) Radial velocity. 

Figure 13. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 4; April 1981. 
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(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 

(b) Radial velocity. 

Figure 14. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 5; April 1981. 
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Nominal altitude = 6.3 km 
F-106B altitude = 6.0 km 





(a) Doppler radar radial velocity contours with F-106B ground track superimposed. Contours are labeled in 
meters per second. 



Elapsed time, min 

(b) Radial velocity. 

Figure 15. F-106B ground track superimposed on Doppler radar radial velocity contours and radial velocity 
at points along flight path. Penetration 6; April 1981. 
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O Penetration 1, figure 10 
QPenetration 2, figure 11 
^Penetration 3, figure 12 
A Penetration 4, figure 13 
(^Penetration 5, figure 14 
CiPenetration 6, figure 15 
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Figure 16. Comparison of Doppler radar and F-106B radial velocity data from April 1981. 
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Radial velocity, m/s 


Altitude, km 


O 3.2 
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Elapsed time, min 
(a) Penetration 3; April 1981. 

Figure 17. Doppler radar radial velocity at three altitudes for three passes through thunderstorm. 


27 










Radar radial velocity, m/s 



Airplane radial velocity, m/s 

Figure 18. Comparison of Doppler radar and F-106B radial velocity data from April 1981 excluding penetrations 
3, 4, and 6. (See fig. 16.) 
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Airplane radial velocity, m/s 

(a) Includes all penetrations for 1980 and 1981. 

Figure 19. Comparison of Doppler radar and F-106B radial velocity data for 1980-1981. 
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(b) Excludes penetrations 3, 4, and 6 of the 1981 data. 
Figure 19. Concluded. 



Appendix 


Determination of Radial Wind Velocity 
From Airplane Measurements 

The wind vector was determined by subtracting 
the F-106B air mass relative velocity from the inertial 
velocity by 

W = V/-V (Al) 

where the inertial components were calculated from 
(fig- Al) 


{ w x t >| (V x .\ ( V cos a cos /3 "j f 0 1 

Wy ' f = \ v y’ r _r \ v sin P cos a > -r < /r > (A2) 
j l V z f ) l V sin a J \ Iq ) 


COS xp COS 0 


T = 


cos 0 sin xp 


L — sin 0 


sin 0 sin <p cos xp 
— sin ip cos <p 

sin ip sin 6 sin <p 
+ cos ip cos <p 

sin <p cos 6 


cos ip cos (p sin 0 " 
+ sin ip sin <p 

sin ip cos 0 sin 0 
— cos ip sin <p 

cos <p cos 0 

(A3) 


The matrix T transforms measurements in the 
F-106B body-axis system to an Earth-fixed or iner- 
tial axis system. (See ref. 10.) The component of W 
in the direction of the Doppler radar is (fig. Al) 


W r = W a (A4) 


lationship. True airspeed is determined from Mach 
number and the acoustic velocity as determined from 
static temperature. (See ref. 11.) 

Laminated balsa flow vanes (described in ref. 12) 
were used for measuring sideslip angle and angle of 
attack. Flight test calibrations were conducted to 
correct the measured angle of attack for upwash ef- 
fects. These data consisted of comparing the mea- 
sured angle of attack with the pitch attitude while 
in a steady, level flight over the airspeed calibration 
range under the assumption of small vertical air mass 
velocities. A second-order least-squares error fit to 
the flight data yielded the following correction for 
angle of attack in degrees: 

a = 0.89506a m - 0.47952 (A5) 

Flow-field corrections for angle of sideslip were not 
used. The vanes were located in the symmetrical lat- 
eral plane and far enough ahead of the airplane that 
sidewash corrections were assumed to be small. For 
these flights, the natural frequency of the vanes was 
estimated at about 40 Hz. The measured natural 
frequency of the nose boom was 20 Hz. To prevent 
aliasing, the a and f3 measurements were filtered be- 
fore they were digitized with a low-pass filter having 
a cutoff frequency of 10 Hz. 

The inertial horizontal velocity components in 
equation (A2) were outputs of the INS. The vertical 
component was determined by integrating the verti- 
cal INS acceleration where 


In equation (A2), airspeed was determined from mea- 
surements of pitot-static pressure and air tempera- 
ture corrected for Mach number effects. The pitot- 
static tube was calibrated by the method in refer- 
ence 11. Position effects were removed based on 
inflight measurements and the ground-camera flyby 
method also described in reference 11. The calibra- 
tion curve for static pressure is shown in figure A2. In 
figure A2 the difference between the measured static 
pressure and the correct static pressure (Ap s ) nor- 
malized by the measured dynamic pressure is plotted 
as a function of measured angle of attack. The cor- 
rected static pressure and measured dynamic pres- 
sure are used in the isentropic perfect- gas-flow equa- 
tion to determine Mach number. The measured total 
temperature, Mach number, and instrument correc- 
tion recovery factor were used to determine static 
temperature from the adiabatic perfect-gas-flow re- 


V z . = -9.807 T A z dt + {V 2 .)o (A6) 

Jo 

The constant (9.807) converts g units to meters per 
second squared. The symbol (V 2 t ) q represents the ini- 
tial vertical velocity, which was determined by equat- 
ing the integrated vertical velocity to the pressure 
altitude change occurring over a time interval T of 
interest, i.e., 


A h p {T) = h p (T) - h p ( 0) = - T V 2 . dt (A7) 

Jo 


and 


W.')o = j; 


-A h p (T) + 9.807 J* J A z dt 


(A8) 
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Figure Al. Body axes (x, y, z) relative to an Earth-fixed system (x', y', z * ) and definition of angle of attack and 
sideslip. The F- 106 B body-axis rates p, g, and r are relative to x 7 , y', and z\ The order of rotation from the 
body to the inertial axis system is The symbol W T represents the component of the wind velocity, 

W, in the direction of the Doppler radar. 
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